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Abstract 

Different mechanisms for the creation of strange non-chaotic dynamics in 
the quasiperiodically forced logistic map are studied. These routes to strange 
nonchaos are characterised through the behavior of the largest nontrivial Lya- 
punov exponent, as well as through the characteristic distributions of finite 
time Lyapunov exponents. Strange nonchaotic attractors can be created at a 
saddle-node bifurcation when the dynamics shows Type-I intermittency; this 
intermittent transition, which is studied in detail, is characterised through 
scaling exponents. Band-merging crises through which dynamics remains non- 
chaotic are also studied, and correspondence is made with analogous behavior 
in the unforced logistic map. Robustness of these phenomena with respect to 
additive noise is investigated. 
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I. INTRODUCTION 



Of the diverse structures that are found in nonlinear dynamical systems, strange non- 
chaotic attractors (SNAs) are among the more exotic. These were first described by Grebogi 
et. al. in a the context of quasiperiodically forced systems, and are characterised by having 
a geometrically strange structure (usually a fractal, which is everywhere single valued but 
piecewise nondifferentiable) with "nonchaotic" dynamics: the largest nontrivial Lyapunov 
exponent A is negative, and nearby orbits do not diverge from each other exponentially. 
Since they were first described a number of characteristics of SNAs have been studied 
theoretically [0-13] and experimentally [[14]- 16]. 



Strange nonchaotic dynamics usually occurs in the vicinity of strange chaotic behavior 
and periodic or quasiperiodic (nonstrange, nonchaotic) behavior. The different mechanisms 
through which SNAs are created, either from regular or from chaotic motion, and the mech- 
anisms through which they disappear, either into regular or chaotic motion — the birth and 
death of strange nonchaotic attractors, so to speak — is a topic of considerable current inter- 
est. 

In this paper, we address this question in the context of a typical dynamical system 
with quasiperiodic forcing and discuss a number of transitions in such systems. These 
include, apart from the transition from quasiperiodic or chaotic motion to strange nonchaotic 
dynamics, transitions between different SNAs. 

We study the quasiperiodically forced logistic map which is given by the equations 

x n+ i = a(l + ecos(27r0„)) x n (1 - x n ) 

4> n+1 = (fi n + lo (mod 1), (1) 

where x E R 1 , (fi E S 1 , u = (\^E — l)/2 is the irrational driving frequency and e represents 
the forcing amplitude. With quasiperiodic driving the periodic attractors of the logistic map 
become quasiperiodic attractors, and following the standard notation || , we denote a torus 
attractor of period n in R^xS 1 as n-T. 



Four mechanisms or scenarios for the creation of SNAs in quasiperiodically driven systems 
have been advanced. In the absence of quasiperiodic forcing, the corresponding systems 
show the by now standard scenarios for the route to chaotic or aperiodic behavior, including 
quasiperiodicity, period doubling and tangent bifurcations, intermittency, crises, and band 



merging or reverse bifurcations |L7 |. There are parallels to several of these in the different 
routes to SNAs. 

• Heagy and Hammel || identified the birth of a SNA with the collision between a 
period-doubled torus and its unstable parent. This mechanism requires that a period 
doubling bifurcation occur, after which the stable torus attractor gets progressively 
more "wrinkled" as the parameters in the system change, i.e. x(4>) becomes more and 
more oscillatory. Following a collision with the unstable parent torus, at an analogue 
of the attractor-merging crisis that occurs in chaotic systems ||, the SNA is born. 
A remains negative throughout the collision process. This mechanism, which we 
denote by HH in the remainder of this paper, has been seen both for a period 2 torus 
as well as for the period 4 torus; the quasiperiodic forcing drives the system into chaos 
well before the infinite sequences of period doubling can occur. 

• The "fractalization" route for the creation of SNAs has been described by Kaneko ||[7| . 
A torus gets increasingly wrinkled and transforms into a SNA without any interaction 
with a nearby unstable periodic orbit. This route to SNA (and eventually to chaos) 



has also been observed in higher dimensional systems |TT 



In the intermittency scenario for the formation of SNAs ||18|| , as a function of driving 
parameter a strange attractor disappears and is eventually replaced by a 1-frequency 
torus through an analogue of the saddle-node bifurcation. In the vicinity of this 
crisis-like phenomenon [[| the attractor is strange and nonchaotic. We have shown 
that the dynamics at this transition shows scaling behavior characteristic of Type I 



intermittency [TI| , and the signature of the transition is an abrupt and characteristic 
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change in the variation of A. The intermittent route is general one which can be found 



in other systems as well [19 



• Yalginkaya and Lai || identify the birth of SNAs with a blowout bifurcation [H] , when 
a torus loses transverse stability. The Lyapunov exponent also has a characteristic 
dependence on parameter in this case. 

SNAs can be quantitatively characterised by a variety of methods, including the estima- 
tion of the Lyapunov exponents and the fractal dimension 0,^,^3, spectral properties 



examination of the time-series [ 1T|J21 1. The geometric strangeness of the attractor can be 



measured through indices such as the phase-sensitivity exponent [10], while the chaoticity 



properties can be studied by examining the finite time Lyapunov exponents fllOylql . 

A number of transitions in this system are investigated wherein three of the above mech- 
anisms for the creation of SNAs are known to be operative (blowout bifurcations cannot 
occur), notably the processes whereby SNAs are created from torus (T) attractors 

...nT«n band SNAs (2) 
2 n T ^ 2 „_i band SNAg ^ _ ^ 

or others, such as from SNA to chaotic attractors (C) or from a fc-band SNAs to a &;/2-band 
SNA, 

. . . n band SNAs <-> n band C .... (4) 
...++k band SNAs jfe/2 band SNAs «-> fc/2 band C++... (5) 

The latter class of bifurcations (for k a power of 2) are studied here through analysis of the 
Lyapunov exponents and their distributions. We also explore the effect of additive noise. 

This paper is organized as follows. In Section II, we describe the "phase-diagram" 
for the forced logistic map: the regions corresponding to the different dynamical behavior 
that obtain are delineated as a function of parameters. This phase diagram is canonical 
for dynamical systems such as Eq. ([!]), namely unimodal maps that are parametrically 



driven, and generalizes the usual bifurcation diagram that is obtained in the absence of 
forcing. Transitions to SNAs and the characteristic behavior of the Lyapunov exponents are 
discussed in Section III, where we also discuss other transformations that are undergone by 
SNAs, such as the analogue of band-merging. The creation of SNAs is often accompanied by 
intermittent dynamics, and this can be quantitatively described in terms of scaling exponents 
at these transitions. These results are also given in Section III, where we discuss, in addition, 
the effects of additive noise. This is followed by a summary in Section IV. 

II. FORCED LOGISTIC MAP: PHASE DIAGRAM 

The quasiperiodically forced logistic map || is particularly convenient for study since the 
phenomenology is smoothly related to that of the logistic map in the limit of e — > 0. Since 
x and are uncoupled in this limit, the period-fc orbits of the logistic map are converted 
into /c-frequency tori, and the chaotic attractors in the logistic map appear as chaotic band 
attractors of the two-dimensional map, Eq. (JJ). 

The region of interest in the phase space is0<x<l,0<</><l. For e 7^ 0, it is clear 
that motion will remain bounded in this region so long as a(l + e cos(27r </>„)) 6 [0,4]. Thus 
for any a < 4, the largest value of e allowed is 4/ a — 1. We therefore redefine the driving 
parameter as e' = e/(4/a — 1) and study the system for < e' < 1. 

Fig. 1 is a phase-diagram |18[ of the system as a function of a and e'. (The regime a < 2 is 
featureless and uninteresting.) The different possible dynamical behavior — periodic, strange 
nonchaotic and chaotic attractors, corresponding to the symbols P, S and C in the phase 
diagram — are characterised through the largest nonzero Lyapunov exponent, 

1 N 

A = lim — Vln|a(l + ecos(27r0 i ))(l -2xi)\. (6) 

AT— >oo N 

i=l 

We calculate A typically from 10 6 iterations of the map (after removing transients for 10 5 
iterations) which is sufficient to converge the results to within 10~ 4 . To obtain Fig. 1, A has 
calculated in a 100 x 100 grid, and the boundary of the chaotic regions are determined from 
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the A = contour. 

The main interesting feature of this phase diagram is the presence of two separate chaotic 
regimes, C\ and C2. The C\ region is the continuation of the chaotic regime in the logistic 
map (this appears, for e = at the end of the period-doubling cascade, at a = 3.5699 . . .). 
The superstable orbit which exists at e = 0, a = 2 continues (slightly distorted) for nonzero 
ef along the locus e' ~ 2(1 — 2/a)/5(4/a — 1); this line (the dotted curve in Fig. la) separates 
Cx and C2. The latter chaotic region is one of low nonlinearity and large-amplitude forcing. 
The two types of chaotic attractors that occur in these two regions are qualitatively quite 
different — see Fig. 2a and 2b. 

SNAs are found in vicinity of the boundaries of the chaotic regions, where A < 0. The 
different mechanisms through which they are formed in different regions are discussed below. 

Along the left edge of the region C2, for 0.59 < e' < 1, the route to SNA is through 
fractalization. This region lies to the left of the curve of "superstable 1-torus" points (the 
dashed line in Fig. la), namely well before torus doubling occurs. For fixed e', SNAs are 
created through fractalization in an interval in a as indicated in the phase diagram. 

Along the right edge of C2, SNAs are born through intermittency. Consider a superstable 
1-torus. In the regime 0.59 < e' < 1, as a is decreased, holding e' fixed, a saddle-node bifur- 
cation takes place and the torus becomes a SNA, eventually becoming a chaotic attractor. 
This process occurs over a somewhat narrow range of parameters as compared to the other 
routes discussed in Section 1. 

Along the left edge of Ci SNAs are created either via the HH mechanism, or through 
fractalization. The phenomenon of truncation or interruption of the period doubling cascade 
resulting from the quasiperiodic driving was noted by Kaneko p2| , who also showed that 
there was a power-law scaling between the number of doublings and the driving amplitude. 
The requirement, in the Heagy-Hammel mechanism, that a period doubled torus collide 
with its unstable parent is a stringent one, and this is achieved only for selected parameter 
intervals. Fractalization is also a likely cause for the interruption of period doubling p2| . The 
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actual mechanism through which a SNA is created is easily identified from the morphology 
of the attractor: a A;-Torus gives rise to a /c/2-band SNA in the Heagy-Hammel mechanism, 
and to a A;-band SNA through fractalization. At low forcing amplitude the fractalization is 
found to be dominant while HH mechanism is more common at higher forcing. In particular, 
we do not find any fractalization event above e' ~ 0.7. 

The regions of SNA are quite complicated. The boundaries between chaotic attractors 
and SNAs are convoluted, with regions of torus attractors interspersed among the regions 
of SNAs. Thus there are several transformations that the attractors of this system undergo 
apart from the birth and death mechanisms, namely the transition from torus to SNA or 
SNA to chaos. Indeed, along a line of constant e or constant e', A is highly nonmonotonic 
and reveals a number of bifurcations. An example is shown in Fig. 3, for constant e = 
0.05. Clearly, there are transitions from T — > SNA — > Chaos — > SNA — > Chaos . . ., finally 
terminating in the Ci chaotic region. 

Within the SNA region, there are transitions from 2™-band SNAs to 2™ _1 -band SNAs. 
At these transitions again there are distinctive and unusual signatures in the dependence of 
A, as well as in the distribution of finite-time Lyapunov exponents. This is discussed in the 
following section. 

Compared to the system at e = 0, for nonzero e it is usually possible to observe regular 
motion only on n-T with n = 1, 2, 3, 4, 6, 8, 16. Higher n-frequency tori are surely present, 
but the regions where they occur are extremely narrow. The period-3 window of the logistic 
map extends as a narrow region of stability wherein one can locate 3-T, 6-T, and the associ- 
ated SNAs, which are formed through the same general mechanisms as outlined in Section I. 
A blow-up of this region of the phase diagram (marked by W in Fig. la) is shown in Fig. lb. 
The intermittent route to SNA that operates on the right edge of C 2 occurs here along the 
top edge of the periodic 3-T region. The 3-band and 6-band SNAs that occur appear to be 
formed through the fractalization or Heagy-Hammel route. 

This qualitative picture appears to be valid for all quasiperiodically forced unimodal 
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maps. We have studied the additively forced logistic map, as well as the forced sine- 
map and found a similar phase diagram. For higher dimensional systems, the picture gets 
somewhat more complicated, but the essential qualitative features carry over, much in the 
same manner as the bifurcation diagram for 1—d maps generalizes for higher dimensional 
dissipative chaotic systems. This system is also not very sensitive to the numerical value of 
the irrational driving frequency u: qualitatively similar phase diagrams have been obtained 
for other cu's. 



III. RESULTS AND DISCUSSION 

In this section we study the variation of the dynamics through several transitions in the 
system as the parameters a and e are varied. In addition to the Lyapunov exponent A itself, 
we also examine the iV-step Lyapunov exponents, A at, namely 

1 N 

X n = T7 E ln + ecos(27r&))(l - 2xi)\, (7) 

iV i=l 

their variance, a, and the distribution 

P(N, X)d\ = Probablity that A at lies between A and A + dX. (8) 

This is of great relevance when studying the stability of systems where a small change in 
control parameter gives raise to drastic change in dynamical behavior. In our calculations, 
A and its variance are typically computed from a sample of 50 estimations of step length 
N = 10 5 . 

The SNA — > chaos route, which has been studied extensively and where it is known that 
the Lyapunov exponent varies linearly through the transition Jl2|, will not be considered 
here. 



A. From Tori to SNAs 
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1. The Heagy-Hammel mechanism 



In the HH mechanism ||, a period 2 n torus attractor gets wrinkled and upon collision 
with the parent unstable period 2 n_1 torus, a 2 n_1 band SNA is formed. This phenomenon 
is similar (in some sense) to the band-merging crisis that occurs in unforced systems, where 
the Lyapunov exponent is known to have a discontinuous slope |p3|. In contrast, in Ref. || 
the transition to SNA appeared to be smooth as a function of the parameter (see Figs. 2 
and 10 in 0). 

We examine this behaviour in some detail in Fig. 4a, where for e' = 0.3 the crisis takes 
place at cyhh = 3.487793 . . .. When examined in a sufficiently small neighborhood of a c , the 
transition is clearly revealed by A: on the torus, A varies smoothly but in the SNA phase 
the variation is rather irregular and the crossover between these two behaviours is abrupt. 
It is also possible to identify the transition point from the examination of the variance in A , 
shown in Fig. 4b: in the torus region, a < olhh-, the fluctuations in A are small, while for 
a > otHH-i the the variance is large and depends irregularly on the function of the control 
parameter. Unlike the case of band-merging or widening crises |23j there is no distinctive 
signature in the variation of the Lyapunov exponent itself (except for a certain irregularity 
in the SNA phase). 



2. The Fractalization mechanism 

During fractalization, a period k torus attractor get wrinkled and eventually forms a k 
band SNA ||12|| . There is no apparent interaction with unstable periodic orbits in contrast 
to the the HH case, and there no analogue of any crisis-like behavior. The variation of the 
A and its variance at such a transition are shown in Figs. 5a and 5b respectively, as a 
changes through t\F = 2.6526 ... at e' = 1. The behavior of variance in Fig. 5b is similar 
to that of the HH case (Fig. 4b), except in the magnitude of the fluctuations. These figures 
also show that the transition from torus to SNA is smooth with no particular signature in 
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A. 



3. The Intermittency route 

Intermittent SNAs are morphologically quite distinct from those formed through other 
mechanisms. An example of the transition to such SNAs is shown in Fig. 6. The behavior 
of A as a is reduced through aj = 3.405808806 ... at e' = 1 (Fig. 6a) is distinctive and may 
be contrasted with Figs. 4a and 5a. The torus which exists for a > aj and the intermittent 
SNA are depicted in Figs. 6b and 6c respectively. On the intermittent SNA, most points 
remain near the parent torus, with sporadic large deviations. On fractalized SNAs points on 
the attractor stay close to the parent fractal torus (Fig. 1 of Ref. 0), while on HH SNAs, 
points are distributed within the entire region enclosed by the wrinkled bounding tori (see 
Fig. 1 of Ref. H). The other characteristic behavior can be extracted from Fig. 6d where we 
plot the variance: this changes abruptly at the transition, where a saddle-node bifurcation 
occurs showing the characteristic signature of the intermittency route to SNAs. 

4- Finite-time Lyapunov exponents 

As is well known, while A is negative on a SNA, for short times, the local Lyapunov 
exponent can be positive. One of the characteristics of the SNAs born through different 
mechanisms is the difference in the distribution of finite-time exponents, namely P(N,X). 
In the limit of large N, it is clear that this distribution will collapse to a delta-function, 
limAT^oo P(N, A) — > S(A — A). The deviations from — and the approach to — the limit can be 
very different for SNAs created through different mechanisms. 

Shown in Figs. 7a- 7c are the distributions for P(50, A) across the three transitions dis- 
cussed above, namely on the tori and corresponding SNAs. A common feature of all three 
cases is that P(N, A) is strongly peaked about A when the attractor is a torus, but on the 
SNA, the distribution picks up a tail which extends into the A > region. This tail directly 
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correlates with the enhanced fluctuation in A on SNAs (see Fig. 4b, 5b or 6d). On the 
fractalized SNA, the distribution shifts continuously to larger A, but the shape remains the 
same for a < aj and a > ai, while on the HH or intermittent SNA, the actual shape of the 
distribution on the torus and the SNA are very different. 

One remarkable feature of intermittent SNAs is that the postive tail in the distribution 
decays very slowly: even for N as large as 10 4 , it does not completely disappear [|TJ|. In 
order to quantify this, we define the fraction of positive local Lyapunov exponents as 



F+(N)= P(N,X)dX, (9) 
Jo 

and similarly, the first moment, namely a local K-S entropy, 

rod 

K(N) = / A P(N,X) dX (10) 
Jo 

Clearly, firnjv->oo F+(N) and liniTv^oo K(N) —>■ 0. Empirically we have found that on 
the intermittent SNA, these quantities show the large N behavior 

F+(N) ~ N~P (11) 

(similarly for K(N)) while for the fractalized or HH SNAs, the approach is exponentially 

fast, 

F + (N) ~exp(- 7 A^). (12) 

The exponents /3 and 7 depend strongly on the parameters a and e'. For the SNAs discussed 
above, we have calculated F + (N) for large A^ (see Fig. 7d) and obtain (3 ~ 0.2, ■jhh ps 0.02 
and 7/ ~ 0.03. K(N) also has a similar slow fall-off at large N for the intermittent SNA 
K(N) ~ N~P' with (3' ~ 0.55, while for the other SNAs, the decay is slow for small A^ but 
exponential for large N. 



B. Scaling at the Intermittency transition 



The intermittency transition from a torus to a SNA is characterised by scaling behavior 
for the dynamics, similar to corresponding behavior in the unforced case . The 'laminar' 
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phase in this case is the torus, while the 'chaotic' phase is the nonchaotic attractor. In order 
to obtain the distribution time in these phases, we co-evolve two trajectories with identical 
(xq, 0o) and e', with different a: since the angular coordinate remains identical, the distance 
between the trajectories is simply the difference in the x n 's. We calculate the time between 
bursts and fit to the scaling form 

r~(a c -a)- e . (13) 

The numerical value obtained for the attractor with a near ai and e' = 1 is 6 — 0.52 ± 0.03 
and e'— 0.65 is 6 = 0.5 ± 0.03 (See Fig. 8a). This suggests that the intermittency is Type-1. 
The Lyapunov exponent itself shows the scaling form [|1^ 



A - A c ~ (a c - a)", (14) 



at fixed e' which can be compared with the probability density [^3, 25, 26 1 in the SNA burst 



phase (Fig. 8(b)). Both these quantities have the same exponent // = 0.37 ± 0.03 at e' = 1. 
Although other SNAs also show intermittency-like dynamical behavior, the scaling form 
Eq. flT4] ) obtains only for the intermittent SNAs, thus providing a means of distinguishing 
these from HH or fractalized SNAs. 



C. Merging of SNAs 

As the parameters a and e are varied, the attractors of the quasiperiodically driven system 
undergo transformation in a manner which is analogus to the undriven chaotic system. 
Similar to reverse bifurcations or band merging in 1— d maps, n-band SNAs transform to 
n/2-band SNAs. Through such a transition, when the dynamics remains nonchaotic and 
strange, A is a good order parameter. Sosnotseva et. al. |fn| , who discovered the first 
example of this transition in the driven Henon and circle maps, demonstrated the merging 
by examining the phase portait. 

Given the fairly narrow range over which SNAs exist in any system, this transition also 
occurs in a restricted range. We find that the SNA bands that are formed by the HH 

12 



mechanism typically do not merge at negative Lyapunov exponent, but only do so at higher 
driving, when they collide in the chaotic region as at a proper analogue of the band merging 
crisis . However, the variation of A at such transitions does not follow a uniform pattern 



as in the unforced case 23,24 



Those SNAs which are formed via fractalization may merge at negative A. Fig. 9a shows 
the variation of A for such an example of a 2-band SNA merging to a 1-band SNA. The 
Lyapunov exponent decreases with increasing nonlinearity. The distribution for short time 
Lyapunov exponents before and after the SNA band merging crisis confirms that the distri- 
bution shifts to lower A; see Fig. 9b. Both these figures indicate that the chaoticity of the 
system (measured either through local (Xn) or global (A) indicators) becomes significantly 
lower after merging crisis. This band merging is different from unforced case |l7|j24l where 
A generally increases as number of bands decreases with increasing nonlinearity. 

D. The Effect of Noise 

An important consideration in the study of SNAs is their robustness. Given the somewhat 
unusual properties of such attractors and the fact that they exist over small regions in 
parameter space, it is natural to examine the effect of fluctuations. This is of particular 
relevance with respect to the experimental observation of SNAs jHirnj . The effect of noise in 



the logistic map has been extensively studied EM , and it is known that noise generally lowers 



the threshold for chaos — systems with additive noise have a larger Lyapunov exponent for 
smaller nonlinearity. Furthermore, transitions and bifurcations get "blurred" in the presence 
of fluctuations. 

To examine some of these effects, we have studied this system through the introduction 
of additive noise in the dynamics, for example as 

%n+i = «(1 + ecos(27T0 n )) x n (1 - x n ) + p£ n (15) 

where p is the noise amplitude and the random variable {£} is 5-correlated in time. As may 



13 



have been anticipated, the addition of noise "smears out" tori, and the threshold values for 
bifurcations typically shift to lower a at fixed e. 

The actual transitions — now from noisy tori to noisy SNAs — survive, and examples of 
this are shown in Fig. 10 for fractalized, HH and intermittent SNAs, where A is shown as 
a function of a with and without noise. Since as a function of decreasing a the intermittent 
SNA is born out of a torus, the effect of noise is to increase the parameter value at which this 
SNA is created (relative to the p = case) while the opposite effect, namely a reduction of 
parameter value for the transition to HH and fractalized SNAs is seen. These behaviors are 
typical: we have verified that SNAs in related systems such as the quasiperiodically forced 
ring and Henon maps behave similarly. 

Upon addition of noise, the global structure of the different types of SNAs remain 
similar — the qualitative behavior such as scaling, band merging etc. are retained — although 
the numerical values of constants and exponents change. The degree of robustness varies 
with nonlinearity, though: for instance, at e' = 1, the intermittent transition survives for 
additive noise of amplitude up to p = 10~ 6 , while at lower e' the transition is robust to even 
larger p « 10~ 4 . The region of scaling in presence of noise is small, and the intermittency 
exponent varies; we have found that the numerical value of 9 (see Eq. (|T3|)) is 0.56 ± 0.04 at 
e' = 1. 



IV. SUMMARY 

In this paper we have described the phenomenology of strange nonchaotic dynamics in 
a prototypical example, namely the quasiperiodically driven logistic map. There are three 
different mechanisms through which SNAs can be created in this system; these routes to 
SNAs have their analogues in the different scenarios for the onset of chaos in dissipative 
dynamical systems |L7 . 



We obtain a "phase-diagram" for the system, delineating the different asymptotic behav- 
ior that are possible as a function of the parameters, namely n-frequency torus attractors, 
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strange chaotic attractors and SNAs. There are two major chaotic regions with different 
characteristics, separated from each other by a region of quasiperiodic and strange nonchaotic 
behavior. As a function of parameters, therefore, the system can show several transitions in 
the dynamics and several of these have been studied. 

To distinguish among the different mechanisms though which SNAs are born, we examine 
not only the manner in which A changes as a function of parameters, but also the variance, 
a. These indices together give a clear indication of the transition from quasiperiodic to 
strange nonchaotic dynamics. We also examine the distribution of local Lyapunov exponents, 
P(N,X) and find that on different SNAs, the fraction of positive exponents F + (N) or its 
moment, the local KS entropy, decay in different manners depending on the type of SNA 
that is formed. 

In addition to three routes that have been decribed previously ]|,|]||, we identify a new 
mechanism for the creation of SNAs. These intermittent SNAs |18| are atypical in that 
they are created at the quasiperiodic analogue of a saddle-node bifurcation fl7| , and the 
signature of the transition is a discontinuous change both in A as well as in a. The chaotic 
component on the intermittent SNA is long lived, giving rise to long chaotic transients: this 
shows up as a slowly decaying positive tail in P(N, A), and a resulting power-law decay for 
F + (N) or K(N). (On other SNAs, in contrast, these quantities decay exponentially.) We 
further characterise intermittent SNAs by establishing the scaling behaviour for residence 
times in the different phases and find that the qualitative picture is in accord with Type-I 
intermittency |T7|]. This behaviour persists in the presence of additive random noise. 

We have described the intermittent SNA in detail and shown that both in its creation as 
well in its morphology, it is distinct from other SNAs, and bears some relation, in terms of 
the phase diagram, to a re-entrant phase. (In the higher dimensional system of the forced 
circle map [T1J for which the phase diagram has been obtained, the intermittent SNA occurs 
in an analogus region). In the vicinity of such attractors, a number of dynamical quantities 
show scaling behaviour. 
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There are other bifurcation phenomena in such systems, and we have examined the 
case of SNA band merging, namely the coalescence of different branches of a multiband 



attractor. At this transition, in contrast to analogous behaviour in the logistic map p3[ , A, 
which remains negative throughout, decreases. A clear understanding of why this is so is 
not available at present. 

In summary, our primary emphasis has been in the use of a number of Lyapunov 
measures — the largest nontrivial Lyapunov exponent, its fluctuations, the distribution of 
finite time Lyapunov exponents and partial moments of this distribution — to characterise 
the different types of transitions to strange nonchaotic attractors that arise in the forced 
logistic map. 

For hyperbolic systems, the theory for the Lyapunov exponent and for Lyapunov mea- 
sures is well developed pTj. For SNAs, the situation is in a less satisfactory state, and our 
work represents initial efforts towards understanding the phenomenology of quasiperiodi- 
cally driven systems. For a number of bifurcations in such systems in general, it is clear that 
the largest Lyapunov exponent is a good order-parameter. It is likely that the considerable 



formalism for such transitions that has been developed for chaotic strange attractors [17.27 



can be applied in large part to the strange nonchaotic regime, but the extent to which the 
theory carries over is an aspect that remains to be explored in future work. 
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Figure Captions 

Fig. 1 (a) Phase diagram for the forced logistic map (schematic). The rescaled parameter e' is 
defined as e' = e/(A/a — 1). T and C correspond to Torus and Chaotic attractors. The 
shaded region along the boundary of T and C corresponds to SNA (marked S). The 
boundaries separating the different regions are convoluted, and regions of SNA and 
Chaotic attractors are interwoven in a complicated manner. The dashed curve marks 
the locus of the "superstable" orbit (see the text). The region of periodic attractors 
can be further demarcated into period 1, 2 and 4 tori as shown (IT, 2T and 4T). W 
denotes the window of periodic behavior corresponding to the period-3 orbit of the 
logistic map. This is shown enlarged in Fig. lb. Intermittent SNAs are found on the 
edge of the C2 region marked I, while the left boundary of C2 has only fractalised 
SNAs. Along the boundary of Ci, both fractalised and HH SNAs can be found, 
(b) A blow-up of the window W indicated in Fig 1(a). This small window shows 
periodic tori of period 3,6,.. and their SNAs that are created either through the 
Heagy-Hammel process or fractalization. As may be expected, other similar windows 
can be seen upon further enlargement. 

Fig. 2 Typical chaotic attractors in the two regions a) Ci, for a = 3.6, e' = 0.5, and b) C 2 at 
a = 3 e' = 1. 

Fig. 3 Variation of A as a function of a for fixed e = 0.05. Note the highly oscillatory 
structure indicative of several transitions in the system. 

Fig. 4 The transition from a period 2-torus to 1-band SNA though the Heagy-Hammel mech- 
anism along the line e' = 0.3 and at oluh — 3.487793 .... a) behavior of A through the 
transition, and b) the variance. 

Fig. 5 The transition from a period 1-torus to 1-band SNA via the fractalization route along 

the line e' = 1 and at a F = 2.6526 a) behavior of A through the transition, and 

b) the variance. 
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Fig. 6 The transition from a period 1-torus to 1-band SNA via the intermittent transition, 

along the line e' = 1 and at aj = 3.405808806 a) behavior of A through the 

transition, b) Plot of the period 1-Torus at a — 3.405809, c) the intermittent SNA at 
a = 3.405808, d) the variance through the transition. 

Fig. 7 Distribution of finite time Lyapunov exponents across the three transitions. Shown 
are the distributions of P(50, A) before and after the transitions, namely on the tori 
and on the SNA. a) Along the HH route, on the torus at a = 3.4874 (□) and on the 
SNA at a = 3.488 (A) for e' = 0.3. b) Along the fractalization route, on the torus at 
a = 2.63 (□) and the SNA at a = 2.66 (A) for e' = 1. c) Along the intermittency 
route, on the torus at a = 3.40581 (□) and on the SNA at a = 3.4058056 (A) for 
e' = 1. d) Variation of F + (N) for the three different SNAs in a), b) and c) showing 
exponential decay in the first two cases, and a power-law decay for the intermittent 
SNA. The respective symbols are D,o and V. 

Fig. 8 Scaling behavior on the intermittent SNA. a) Plot of the average time between bursts 
vs (a c — a) at e' = l(o) and at e' = 0.65(D). (b) The probability density (j>b) in the 
burst phase (o) and A (□) vs (a c — a) at e' = 1. The measured exponents are // ~ 0.37. 

Fig. 9 a) Variation of A at band-merging bifurcation in SNA for 2-band SNAs to a 1-band 
SNA along e = 0.05. b) Distribution of local Lyapunov exponents, P(100, A), across 
the merging transition in a) which takes place at a c xs 3.387439. 

Fig. 10 The shift in A with additive noise near the threshold for the transition from SNA to 
chaos for the cases of a) HH, b) Fractalization, for p = (solid line), 10~ 3 (dashed 
line) and c) Intermittent SNAs, p = (solid line) p = 10~ 5 (dashed). All quantitites 
are estimated from a set of 50 samples of 10 5 time-steps. 
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